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Abstract: We describe a method for evaluating chiral gauge theories that is not plagued by the doubling problem. To 
demonstrate the efficiency of the approach, we apply our ideas to the chiral Schwinger model. 



1. INTRODUCTION 

Chiral gauge theories play an important role 
in particle physics. Yet very little is known about 
their properties. The main reason for this is the 
lack of a calculational method which takes into 
account non-perturbative effects. 

The lattice regularization was invented for 
that purpose. Until now we have, however, not 
been able to find a consistent lattice formula- 
tion of these theories. Due to the doubling prob- 
lem one is always led to an equal number of 
left- and right-handed fermions with net chi- 
rality zero. Attempts to remove the unwanted 
"doublers" by pushing their masses to the cut- 
off while preserving chiral symmetry have, in 

* TALK GIVEN AT THE WORKSHOP ON NON- 
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spite of great efforts, not brought the desired 
success [Q. Some authors have arranged them- 
selves with this situation by postulating the exis- 
tence of mirror fermions ||^ . But even these mod- 
els are inadequate. A non-vector theory leads - 
after the anomalies have been cancelled - to a 
gauge invariant, complex effective action whose 
imaginary part is non-zero Q . The action of the 
mirror-fermion model is, however, real. 

But it may also be that non- vector theories do 
not exist as consistent field theories with non- 
trivial dynamics. A first answer to this ques- 
tion can be obtained from the chiral Schwinger 
model. In four dimensions one has to resort to 
numerical methods to address this problem. We 
believe, however, that it should at least be pos- 
sible to answer some of the more fundamental 
questions, in particular whether the imaginary 
part of the effective action is indeed non-zero. 
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In order to evade the doubling problem, we 
propose to evaluate the fermionic action in the 
background of continuum gauge fields, which we 
will construct fi'om lattice gauge field configu- 
rations. The starting point is a formula for the 
effective action, which was derived some time ago 
by Alvarez-Gaume et al. for a similar purpose. 
The virtue of this formula is that it involves the 
ordinary covariant Dirac operator only and that 
it is amenable to explicit calculations. The basic 
idea of the method as well as some theoretical 
background are presented in sec. 2. In sec. 3 we 
sketch the derivation of continuum gauge fields 
from lattice gauge field configurations. Our ideas 
are then applied to the chiral Schwinger model 
in sec. 4. Finally, in sec. 5 we present some con- 
cluding remarks. 



2. EFFECTIVE ACTION 

We will work in Euclidean space compactified 
to the torus T^. But one can also think of other 
geometries. Consider now the Dirac operator 



1 T75 



(1) 



which minimally couples left-handed and right- 
handed fermions, respectively, to the gauge fields. 
The gauge fields are assumed to be in a complex 
representation of some gauge group G, 



A — AO-rpa /IT — _ 4 



(2) 



We are interested in the effective action for left- 
handed fermions, S-.{A). This is given by 



(3) 



It is understood that eq. (Q) and the following ex- 
pressions are properly regularized, without com- 



mitting ourselves to any particular scheme for 
the moment. For the real part of S-{A) one finds 

2ReS'_(A) = S-{A) + S+{A) ^ S{A) + ^(O) (4) 

with 

e-^(-^) = j Vi^Vije-^'^^^^^, (5) 

(6) 



tp{A)=t{^+4), 



which means that it is given by the action of 
the corresponding vector theory. This part of the 
action does not constitute any problems. It can 
be computed by standard lattice techniques. 

Thus the chiral nature of the fermions is re- 
flected entirely in the imaginary part of the ef- 
fective action. Following Alvarez-Gaume et al. 
we shall now derive an expression for ImS*- (A) 
on which we will base our further discussion. Our 
derivations will be purely formal. For rigorous 
proofs the reader is refered to the original liter- 
ature p|. We start from the expression 



S.{A)-S.iO) 

= -In det i ^_ (A) + In dot i ^_ (0) 
= -Tr In i ^_ (A) + Tr In i (0) 
= -Tr In (1 + {i ^)-^ijfiP-) . 
Expanding the logarithm, this gives 
5_(A)-5_(0) 



(7) 



(8) 



By replacing Ahy^A — tA on the right-hand side 
of eq. (^), differentiating the whole expression 
with respect to t and integrating over t again, 
one obtains after doing the sum over n again 
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S-{A)-S-{0) 



f 

Jo 



(9) 



dtTTi4{i^{*A)y'p^ 



This results in the expression for the imaginary 
part 

ImS-iA) = -ll'dt Tr 75 PMr^ (10) 



While the real part of the action is gauge in- 
variant, the imaginary part is generally not. Un- 
der gauge transformations 



it transforms like 

lmS-{<^A) - lmS-{A) = 2ttQ5, 



(11) 



(12) 



where Q5 is the Wess-Zumino action. If the the- 
ory is anomaly-free, Q5 — and S^{A) is gauge 
invariant. But we still expect lmS^{A) ^ 0, also 
in a model with mirror fermions. 

A straightforward lattice calculation of eq. 
( p^ ) with, e.g., Wilson fermions would give a di- 
vergent result in the continuum limit. The sit- 
uation is similar to the case of the topological 
charge. 



Q 



167r2 



d^x Tr FF, 



FF — p(T F F ^ 



(13) 



(14) 



where the naive latice transcription of FF |4j 
does not have any topological significance. In- 
deed, it has been shown |^ that the gauge invari- 
ant part of ImS'_ {A) is given by the ry- invariant 
of the five-dimensional Dirac operator. 



which is a topological invariant. 

Let us stay with the example of the topological 
charge for a moment. It is well known that the 
naive expression for Q receives perturbative con- 
tributions, which are divergent in the continuum 
limit and so cover all topological effects. The ori- 
gin of these contributions lies in fluctuations of 
the gauge fields on the scale of the lattice spac- 
ing. These fluctuations are unphysical near the 
continuum limit, and only those theories which 
have a continuum limit are considered here. If we 
now extrapolate the lattice gauge field smoothly 
to the interior of the hypercubes - thereby re- 
specting the fact that a topologically non-trivial 
gauge field configuration cannot be in the same 
gauge everywhere on the torus - and calculate Q 
in the background of this continuum gauge field 
using the expression (p^, we find not only that 
all divergences have been eliminated ^ , but also 



H = i-i^dt+ipCA), 



(15) 



that Q is an integer as it should. (See the next 
section for further details.) This is the essential 
idea behind the geometrical expressions for the 
topological charge |^J^,^ , though the actual cal- 
culation of Q proceeds via the transition func- 
tions, which connect the gauges in the different 
hypercubes, and for those it is sufficient to know 
the gauge field on the boundary of the hypercubes 
only. Given the continuum gauge field, it would 
equally well be possible to compute the topologi- 
cal charge by discretizing this field again, but on 
a finer mesh, computing the corresponding link 
matrices and using the naive lattice expression 
for FF Q to obtain Q. The continuum value of 
the charge is then found by making the mesh 
finer and finer until the result has converged. 

We now propose to compute Im5'_(A) in an 
analogous way. We start from a lattice gauge 
field configuration, which we regard for the mo- 

We assume that the action has been chosen such that 
the susceptibihty is not affected by dislocations H] . 
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ment as being given to us. The first problem is 
then to extrapolate the lattice gauge field to the 
whole interior of the hypercube. This construc- 
tion, which is highly non-trivial, will be sketched 
in the next section. Further details can be found 
in the original publication The next prob- 
lem is to compute luiS-{A) in the background 
of this field. As we have already indicated, the 
expressions for the effective action have to be 
properly regularized in order to make sense. We 
suggest to use the lattice regularization with a 
lattice spacing, a, much smaller than the origi- 
nal lattice spacing. The latter we shall set equal 
to one for convenience. For Wilson fermions we 
then find 

ImS^iA) ^\fdt a4^{Tr 75(7^. - 1) 

X ^In U^^,, Ul^g^+f,.^ + Tr 75(7^ -f 1) (16) 

where U are the link matrices on the fine sublat- 
tice, and Q is the inverse fermion matrix, 

-^E(^A< + 1)^."-A,M^-A,. (17) 

4 1 

a a* 

It is straightforward to compute the link matri- 
ces Z//, as well as the logarithm and powers of U, 
from the continuum gauge field. The right-hand 
side of eq. (^6|) can be very efficiently calculated 
by means of a stochastic estimator . This re- 
quires only one inversion of the fermion matrix 
for one set of random numbers. In the final step 
we then have to perform the calculations on finer 
and finer sublattices until we obtain a stable re- 
sult (if this exists). 



We like to emphasize that eq. (^6|) is only one 
out of many possibilities of computing Im S- {A) 
in the background of a continuum field. An al- 
ternative method would be to use the gauge field 
A^{x) directly instead of In Ux,^i. U^.^^. But this 
would require a further regularization step. One 
such step is the point-splitting technique. Which 
method is best suited can only be decided in con- 
crete calculations which we will return to in the 
future. 

It may be possible that there are more elegant 
ways of computing the imaginary part of the 
effective action. Alvarez-Gaume has attempted 
a geometric construction similar to that of the 
topological charge, but without success [0. In 
two dimensions it turns out that one can calcu- 
late the effective action analytically, as we will 
show in sec. 4. 

Whether we can do simulations of chiral gauge 
theories depends now on the magnitude of Im5_ {A) . 
So far we do not know how to handle complex 
actions. If ImS*- {A) is small, we may simulate 
gauge field configurations according to the real 
part of the effective action (including the gauge 
field action) on the lattice in the standard fash- 
ion and include the phase factor in the statistical 
average over the ensemble. 



3. CONTINUUM GAUGE FIELDS 

We shall now sketch the construction of con- 
tinuum gauge fields Af^{x) from lattice gauge 
field configurations. We work on a hypercubic 
lattice with periodic boundary conditions. The 
lattice points are denoted by s, and the link ma- 
trices are called Us,ii- 

The gauge fields have to fulfill the following 
constraints, (i) The link matrices derived from 
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Af^ must agree with those of the original lattice, 
Us^fj,. (ii) The construction must be gauge covari- 
ant, which means that a lattice gauge field trans- 
formation must result in a continuum gauge field 
transformation of A^{x). (iii) The gauge field 
Afj_{x) should lead to the fiber bundle which has 
been previously derived from the lattice gauge 
field configuration This guarantees that the 
topological charge is the same in both cases. 

In the first step we construct in each hyper- 
cube, 

c(s) = {x e I < x^ <s^ + (18) 

a local gauge field A^ff\x), such that on the in- 
tersection of two adjacent hypercubes the corre- 
sponding gauge fields are connected by a gauge 
transformation given by Liischer's transition func- 
tions p . In the second step we then construct a 
global, and generally singular, gauge field out of 
the local gauge fields. 

Because the construction involves derivatives 
of the transition functions in intermediate steps, 
it is necessary to enlarge the hypercubes to open 

sets, 

c(s) = {x e T"* I Sf,-e <x^< Sf, + l+e,yfi}Xl9) 

and define smeared transition functions on open 
faces /(s, /i) = c(s) n c(s — fi): 



Vs,tJ.ix) = Vs,^{si + (j){xi - Si), ■ 
• • • , S^, • • • , S4 + - S4)). 



(20) 



Here Vs^f_i is Liischer's transition function [|6| de- 
fined on /(s, /i) = c{s)r\c{s—fl) and is a smooth 
function, : R ^ [0, 1], with the properties 



for t < 77 

1 for t > 1 — T] 



(21) 



and e < 77 < 1/2. It follows that the w's satisfy 
the cocycle condition 



Vs-il,,^{x)Vs,f^{x) = Vs-0,fj.{x)Vs,i'{x) 



(22) 



on p{s,^,v) = f{s,^) n f{s — v,^) and there- 
fore define a principal bundle over T"* . The gauge 
field a\!\x) can now be constructed step by step 
making use of the relation 



A'-r'\x)^isAx)A'^\x)i^Ux) 



+Vs,iy{x)di_,v^l{x). 



(23) 



Starting from the value zero at x G c(s), — e < 



x„ ~ < e for Vz/, eq. (|23|) gives the gauge field 
in the neighborhoods of the other corners of c(s), 
where it is also zero. In the next step one takes 
the gauge field to be zero in e-neighborhoods of 
the links in c(s) which originate from s and ap- 
plies eq. ( p^ ) again. Then one defines A|f (x) on 
the six plaquettes p(s, /x, v) by interpolating the 
expressions already obtained. By means of eq. 
( p3| ) one finds A\f\x) on all the other plaquettes 
contained in c(s) . Continuing in this fashion one 
arrives at the expression 

^''^Hx) = (piXa - Sa)<i){xp - Slj)(j){x^ - S^) 
X {Z^{s, a, /3, 7 I Sa -I- 1, S/3 + I, + 1, X^) 



XWs+a+;§,;3)(S" + 1: S/3 + 1, X^, Xf,)) 

XZ^{S + a + P,"/ \ Sa + 1,3/3 + 1,3^ + I, Xf^) 

+Ad{Vs+a,a{Sa + 1, X0, X^, X^) 

X^s+q+/3./3('Sq + f : S/3 + I7 2^71 Xi_i)) 

XZ^{S + a + (3,"f \ 3a + l,3[3 + l,3y + I, X ^) 
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^Vs+a+j./y{Sa + Ij Xp, + 1, X^)) 

xZfj,{s + a + ^, P \ Sa + 1,3/3 + l,Sy + l,Xf,) 

-Ad(-D5+Q,a(sQ + l,Xi3,X^,Xf_,)) 

xZ^{s + a,/3,7 I Sa + 1, s/3 + 1, + 1,2;^)]} 

+ X] - Sa)(j){xp - Sp) 

cycl . perm . 

xZ^{s,a,P I Sa + 1, S/3 + 1,X-^,X^) 
+ (f>{Xa - Sa)4>{xp - Sfj) 
xAd{Vs+a,a{Sa + 1, X p , X^ , X f,)) 

xZf^{s + a,(3 \ Sa + l,sp + l,Xy,Xfj^) (^^) 

+ (/)(Xq — Sa)(f>{x^ — S^) 
xAd{Vs+&^aiSa + l,Xp,X^,Xf_,)) 

xZf,{s + a,7 I s„ + l,xp,Sj + 
+(j){xa - Sa)Zfj.{s, a\ Sa + l, xp, x^, a:^)] , 
where 
Z^(s,a,/3,7 I x) 



due to the cocycle condition. It is straightforward 
to verify that for x G /(s,/i) 

P exp{ /' dtA^^^ {x+{l- = (26) 



where ^_|_^ are the Unk variables in the com- 
plete axial gauge defined in ref. ||]. In eq. ( |2^ ) 
one can take the limit e — > and put (f)(t) — t 
for < < < 1. 

In order to obtain the global gauge field we 
make use of the section , originally defined on 
the boundary of the hypercube JT^ and extended 
to its interior in ref. H, and write 



A^,{x)^w'{x)-^a'^^\x)w'{x) 



(27) 



+w^{x) ^d^w^{x). 

Under a lattice gauge transformation, 

U{s, ^ U{s, m) = g{s)U{s, y)g(s + fir\ (28) 

the local gauge field transforms as 

A''^\x)^g{s)-'A^;;Hx)g{s), (29) 

and the section transforms as 

w'{x)^g{s)w'{x)g{xr\ (30) 

where g{x) is the lattice gauge transformation 
extrapolated to the interior of all the hypercubes. 
The latter is given in ref. This gives the de- 
sired result 



Zf,{s,a,(i\x) -{i,+&^^v^^^^pfj (25) A^{x) ^ g{x){A^,{x) + d^,)g{x)- 

Z^(s,a I x) = {is+&,ad^,v~lg^ J{x), 
kd{v)M = iMv-^. 



(31) 



Note that Z^(s,q;,/?, 7 | x) and Zfj_{s,a,P \ x) 
are symmetric in a,/?, 7 and a,P, respectively, 



4. AN EXAMPLE: THE CHIRAL 
SCHWINGER MODEL 

We shall now apply our method to the chi- 
ral Schwinger model. A lot is known about this 
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model from other studies . Our approach dif- 
fers from most of the other approaches by work- 
ing in Euchdean space. Questions regarding the 
positivity of the transfer matrix and the hke will 
and can be ignored here. 

The action of this model is 



Our choice of the 7 matrices is 
/ 1 \ ( ~i 



(32) 



1 



71 



with 



75 = «7i72 = 



,72 



-1 



(33) 



(34) 



The gauge fields will be written = iA^^, where 
A° is real. This then leads to the Dirac operator 



i ^{A) = i 



-52 - A? 



di+A°2 



92 -A? 



(35) 



We may assume that i p{A) has no zero modes 
due to the vanishing of det i p{A) in this case. 
Being a two-dimensional model, it is easy to find 
a closed form for the inverse Dirac operator in 
eqs. (p|jlO|), which is the solution of 



z piA)iz p{A))-^x,y) = 5\x - y) 1. (36) 
We obtain 



{zP{A))-\x,y) = 
.4 



(27r)2 



(37) 



where 

/(x, y) = exp{/ (Pz{s{x, z) - s{y, z)) 
x(A?(z) + zA§(z))} 



(38) 



with 



s(x, z) 



Jq{x-z) 



-(-(71 +iq2). 



(39) 



This result generalizes to other gauge groups as 
well. 

We are primarily interested in the gauge in- 
variant sector of the model, as it arises e.g. af- 
ter the anomalies have been cancelled. The (non- 
trivial) dynamics in this sector is driven by the -q- 
invariant. To keep the calculations in this write- 
up short, we shall fix the gauge and state the 
general result only at the end. We choose Lan- 
dau gauge, d ■ A^{x) = 0. In this case eq. (|3S 
simplifies to 



/(x,y) = exp{i J — 
x{~q,A^,{q)+q,AM)}. 







(40) 



On the torus (or any other compact space) the 
momenta are discrete. So, whenever we write an 
integral over internal momenta, the appropriate 
sum respecting the boundary conditions should 
actually be understood. In our case it is, how- 
ever, simpler to do the integrals rather than the 
corresponding sums. In the infinite volume limit 
both methods should give the same result. 

We are now ready to compute the effective ac- 
tion. We choose the point-splitting method to 
regularize the integrals, which is computation- 
ally simple. Using the Pauli-Villars regulariza- 
tion gives the same result. It would have been 
even more appropriate to use the lattice regu- 
larization as proposed in sec. 3, in particular as 
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the continuum gauge fields can be stated explic- 
itly in this case However, it is hopeless to 
find a closed form for the inverse of the lattice 
Dirac operator, which means that one has to do 
the sums numerically. We shall return to that - 
as a check of the method - when we will apply 
our ideas to four-dimensional models. But we be- 
lieve that the results will be the same. Writing 
now *A = itA'^, we have 

ImS-(A)^ lim -Xt f dt f (fx 



xTrz75 40(a;) 



(27r)2 fc2 



(41) 



where 

(I 



£2 



(42) 



The integral over k can be done immediately and 
gives 



27r {x — yy 



(43) 



For the integral over y only the region y x 
is relevant. We may therefore expand the line 
integral and f{x,yY about y = x. This gives 



exp{i /J" ds^^tAl) = 1 -f i{xi - yi)tA\{x) 
+i{x2 - y2)tAl{x) + ■ ■ ■ 
and 

f{x,yy ^\ + {xi-yi)tAl{x) 
-{x2 ~ y2)tA\{x) + ■ ■ ■ . 



(44) 



(45) 



If we insert this now into cq. (ETl) we obtain 



lmS-{A) = — 3- / dtt I d^xTv ij5 
47r Jq 



A0(a;)2 



(46) 



A°{xy 



= 



(A"(x)2 = A'({x)^+A^{x)^). It is straightforward 
to repeat the calculations in an arbitrary gauge. 
In momentum space we find 



Im._(A) = ^/ 



1 f d^q 1 



(27r)2 g2 

x{q,AUq)~q2AUq))q-A^{q), 



(47) 



which is the anomalous term. Altogether this 
means that the gauge invariant sector of the 
chiral Schwinger model is identical with the 
corresponding vector theory, i.e. the massless 
Schwinger model. In particular it says that the 
ry-invariant is zero. 

Along the same lines we can compute the real 
part of the effective action, 

Re{S^{A) - 5_(0)) = i(5(A) - 5(0)). (48) 

This amounts to replacing 75 by il in eq. (^6|), 
which gives 



S{A) ~ 5(0) = dtt j d^ 



xTr 



A^{x) 



A°{x) 



(49) 



^ j d^xA\xf 



This is the well known result, which says that 
the gauge field acquires a mass 

= — . (50) 
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Note that in our notation the pure gauge field 
action reads ^ / cPxF^^^. It assures us that our 
calculations are correct. In the chiral Schwinger 
model the mass is = e^/27r because of eq. 
(p8|). In an arbitrary gauge eq. ( p9[ ) is found to 
be supplemented by a pure gauge term. 
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5. CONCLUSIONS 

We have formulated a method, which should 
allow us to answer some of the principle ques- 
tions in chiral gauge theories after years of fruit- 
less attempts. 

We have confronted our ideas with the chi- 
ral Schwinger model and found that the gauge 
invariant sector of this model equals the vector 
Schwinger model. This basically means that the 
chiral Schwinger model has no non-trivial dy- 
namics. 

In four dimensions the main task is now to 
evaluate eq. ( p^ ) (or any other equivalent ex- 
pression for lmS-{A)) and demonstrate that the 
result converges as a ^ 0. In a model which is 
anomaly-free we expect the result to be gauge 
invariant. If the result for ImS'_(^) turns out to 
be small, we can then think of numerical simula- 
tions. If, however, lmS-{A) is large and strongly 
fluctuating, we do not see any chance for lattice 
simulations at the moment. We like to emphasize 
that the problem of complex action is inherent 
in any other approach. Whether we can make a 
statement as to whether ImS"- (A) is non-zero or 
not remains to be seen. 

Finally, we like to say that it might still be 
possible to find a "geometrical" expression for 
luiS-{A) or the ry-invariant, and one should 
search for it. This would improve the situation 
considerably. 



References 

[1] For a recent review see, e.g., 

M. F. L. Golterman, Nucl. Phys. B (Proc. Suppl.) 
20 (1991) 526; 

see also the talks at this Workshop 
[2] I. Montvay, Phys. Lett. 199B (1987) 89; 

L. Lin, L Montvay, H. Wittig and G. Miinster, Nucl. 

Phys. B (Proc. Suppl.) 20 (1991) 601; 

see also the talk by I. Montvay at this Workshop 
[3] L. Alvarcz-Gaume, S. Delia Pietra and V. Delia 

Pietra, Phys. Lett. 166B (1986) 177; 

L. Alvarez-Gaume and S. Delia Pietra, in Recent 

Developments in Quantum Field Theory, eds. J. 

Ambjorn, B. J. Durhuus and J. L. Petersen, p. 95 

(North-Holland, 1985); 

S. Delia Pietra, V. Delia Pietra and L. Alvarez- 
Gaume, Commun. Math. Phys. 109 (1987) 691 

[4] P. Di Vecchia, K. Fabricius, G. C. Rossi and G. 
Veneziano, Nucl. Phys. B192 (1981) 392; 
K. Ishikawa, G. Schicrholz, H. Schneider and M. Te- 
per, Phys. Lett. 128B (1983) 309 

[5] M. Gockeler, A. S. Kronfeld, M. L. Laursen, G. 
Schicrholz and U.-J. Wiese, Phys. Lett. 233B 
(1989) 192 

[6] M. Liischer, Commun. Math. Phys. 85 (1982) 29 
[7] A. PhiUips and D. Stone, Commun. Math. Phys. 

103 (1986) 599 
[8] M. Gockeler, M. L. Laursen, G. Schicrholz and U.-J. 

Wiese, Commun. Math. Phys. 107 (1986) 467 
[9] M. Gockeler, A. S. Kronfeld, G. Schicrholz and U.-J. 

Wiese, HLRZ preprint 92-34 (1992) 
[10] See, e.g., 

K. Bitar, A. D. Kennedy, R. Horsley, S. Meyer and 

P. Rossi, Nucl. Phys. B313 (1989) 348; 

M. Gockeler, R. Horsley, P. Rakow, G. Schicrholz 

and R. Sommer, Nucl. Phys. B371 (1992) 713 
[11] L. Alvarcz-Gaumo, private communication 



10 



M. Gdckeler and G. Schierholz / Chiral Fermions 



[12] M. Gockeler, A. S. Kronfeld, M. L. Laursen, G. 

Schierholz and U.-J. Wiese, Nucl. Phys. B292 

(1987) 349 
[13] see, e.g., 

R. Jackiw and R. Rajaraman, Phys. Rev. Lett. 54 

(1985) 1219; 

H. O. Girotti, K. D. Rothe and H. J. Rothe, Phys. 
Rev. D33 (1986) 514; 

N. K. Falck and G. Kramer, Ann. Phys. 176 (1987) 
330; 

D. Boyanovsky, Nucl. Phys. B294 (1987) 223; 

C. Aragao de Carvalho, K. D. Rothe, C. A. Linhares 

and H. J. Rothe, Phys. Lett. 194B (1987) 539; 

T. Berger, Cambridge preprint DAMTP-92-02 

(1992) 

[14] J. Schwinger, Phys. Rev. 128 (1962) 2425 



